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The global geometric entanglement is studied in the context of newly-developed tensor network algorithms 
for finite systems. For one-dimensional quantum spin systems it is found that, at criticality, the leading finite- 
size correction to the global geometric entanglement per site behaves as b/n, where n is the size of the system 
and b a given coefficient. Our conclusion is based on the computation of the geometric entanglement per spin 
for the quantum Ising model in a transverse magnetic field and for the spin-1/2 XXZ model. We also discuss 
the possibility of coefficient b being universal. 
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Introduction. At zero temperature, quantum many-body 
systems exhibit important collective phenomena. The com- 
plex structure of these systems has attracted great attention 
in recent years. In particular, the role played by quantum 
correlations (or entanglement) in quantum phase transitions 
has been the subject of numerous studies |l{]. Several re- 
sults for the ground states of quantum many-body systems 
in one spatial dimension (ID) show that a number of entan- 
glement measures obey universal scaling laws close to or at 
quantum critical points |2j-|6|]. These results are, in part, due 
to the unique properties of the reduced density matrices of 
the ground states ||6[] . Complementarily, the so-called fidelity 
approach to quantum phase transitions BT HlOll has proven spe- 
cially fruitful in determining the existence and nature of pos- 
sible critical points in quantum many-body systems. 

There has also been a growing interest in investigating the 
distance between the ground state of a quantum many-body 
system and the closest separable state in the Hilbert space. 
This idea has been quantified in terms of the so-called ge- 
ometric entanglement (GE) fill which is a measure of 
the global quantum correlations of a quantum many -body sys- 
tem. Remarkably, the study of the GE in ID many -body sys- 
tems has proven fruitful because of its connections to both 
the Renormalization Group and Conformal Field Theory [5]. 
Moreover, the study of the GE allows to identify separable 
ground states of quantum many -body systems 11311 . As em- 
phasized in Refs. lfl3i[l4ll . the occurrence of such factorized 
ground states is a typical signature of quantum phases of mat- 
ter with symmetry-breaking order, as opposed to phases with 
topological order where the ground state is globally entangled. 
While these properties show the conceptual significance of the 
GE, the GE has unfortunately also proven to be very hard to 
compute both analytically and numerically. In the case of ID 
systems, many previous studies of the GE have focused on its 
properties for infinite size systems, where translational invari- 
ance together with the infinite-size limit can be fully exploited 
to simplify the calculations. Its behavior for finite size ID sys- 
tems has only been analyzed in a few special cases 1U2I1 . 

In this paper, we show that the GE can be easily extracted 
for finite systems with the aid of newly developed Tensor 
Network (TN) algorithms to simulate strongly-correlated sys- 



tems. This connection allows us to further extend the study of 
the GE by obtaining the finite-size corrections to the GE for 
two important ID models at criticality: the quantum Ising and 
XXZ spin-1/2 chains. The calculations consist of two steps: 
first, we compute a Matrix Product State (MPS) representa- 
tion of the ground state using TN algorithms. Second, we de- 
termine the GE from the MPS representations thus obtained. 
Our numerical observations can be summarized as follows: 
for a ID system at criticality, the finite-size corrections to the 
GE per site £„ for a system of size n obey 
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n » 1, 
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with b a given coefficient. As we shall see, &„ can also be 
interpreted as twice the free energy per site of a classical ID 
vertex model. The possibility of coefficient b being universal 
is also briefly discussed. 

Global geometric entanglement per party. Let us now in- 
troduce the measure of entanglement that we shall use in this 
paper. Consider a pure quantum state of n parties G H — 
(££) /=1 < H^ > where is the local Hilbert space of party i. 
The global multipartite entanglement of can be quantified 
by considering the maximum fidelity A max between the quan- 
tum state \ifr) and all the possible separable and normalized 
states \<p) of the n parties 



A max = max \(if/\(p)\. 



(2) 



A well-defined global measure of entanglement is obtained by 
taking the base-2 logarithm of A 2 nax , 



£0A) = -io g2 A;L 



(3) 



Here, we will be interested in the above quantity per party, 
which is given by 



& n - n E(i[f). 



(4) 



The above quantity constitutes what we call the global GE per 
party, and has been investigated in different contexts Jf, U, 
12L LL6D . Here we study the finite-size properties of S„ as the 
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size n of the system changes. Finally, we choose each party to 
be a single spin, and therefore & n corresponds to the GE per 
spin. 

Geometric entanglement from tensor network algorithms 
for ID systems. For ID quantum many-body systems, sev- 
eral methods have been developed to compute their ground 
state properties based on MPS representations of their wave 
functions. Here, we use two of these algorithms to compute 
ground states of ID quantum many -body systems under Pe- 
riodic Boundary Conditions (PBCs). These algorithms are 
based on the variational optimization of the expectation value 
of the energy using (i) a direct variational method B17I1 . and (ii) 
a variational Monte Carlo method [18]. We also show how 
to numerically extract the GE once the MPS representation 
of the ground state wave function has been obtained. In all 
our derivations we assume that the Hamiltonian for n sites is 
made of nearest-neighbor interactions, H = YII=\ h^ 1, \ with 
n + 1 = 1 under PBCs and h^ h • ^ a nearest-neighbor two-body 
Hamiltonian. 

Let us first recall some of the key ingredients of the direct 
variational algorithm from Ref. 111711 to find ground states of 
ID many-body systems with PBCs. We start from an initial 
guess state \tffo), expressed in terms of an MPS as 

l<Ao>= J] Tr[A [11 (^i)A [21 ( i2 )...A [ " ] ( i „)]| il ,--- ,s„), (5) 

where s,- = 1 , . . . , d for i = 1 . . . , n, and A w (s;) are D x D ma- 
trices, with d denoting the dimension of the Hilbert space at 
every site i and D denoting the so-called bond dimension. The 
variational optimization algorithm finds an MPS approxima- 
tion of the ground state of the system by iteratively optimiz- 
ing the expectation value of the Hamiltonian H. Specifically, 
at every step the algorithm finds the optimal MPS matrix at a 
given site. For instance, focusing on site m, the expectation 
value of the Hamiltonian can be written as 

E = (mm = Alm|t ^ Alml (6) 

<<AI<A> /Ctat N m „ A>1 ' 

eft 

where vector A^"' 1 has the components of matrices A [m, (s m ) 
arranged as a vector, and and N'^ are effective matrices 
that depend only on the Hamiltonian H and the MPS matrices 
A w (s,) for i + m. The update of the d matrices A w (s m ) at site 
m follows from minimizing Eq. (O with respect to A^ m] . This 
amounts to solving the generalized eigenvalue problem 

fl« A lm] = viV™ AH (7) 

with v the expectation value of the energy. By updating the 
MPS matrices site by site, the approximate ground state is 
obtained as the energy E monotonically decreases. 

Next, we review the basics of the variational quantum 
Monte Carlo optimization algorithm (VQMC) from Ref. IU8I1 . 
This algorithm takes advantage of translational invariance as 
well as of other symmetries in the system to accelerate con- 
vergence. More specifically, for a translationally invariant 



system, the MPS of its ground state is described by just a 
single site-independent set of d matrices, that is, A [!] (s,) = 
A(s) Vi. Furthermore, for systems invariant under reflection 
symmetry, the condition A T = A is imposed. The algo- 
rithm starts by giving random values to the components of 
matrices A (s). Then, these components are changed in or- 
der to minimize the ground state energy. The changes in the 
components A(s) a p at step k in the algorithm follow the rule: 
A(s) aP -> A(s) afi -6(k)- r(s) afi -sign(dE /dA(s) a/3 ), where r(s) afj 
is some random number between and 1 , and 6{k) is the maxi- 
mum allowed variation at step k (which is set to monotonically 
decrease with k). Thus, each component A(s) a p is modified 
independently by a random and well bounded variation, and 
the acceptance or rejection of this variation follows a usual 
Monte Carlo rule. Also, as the expectation value of the en- 
ergy decreases and converges, more sampling is required to 
avoid undesirable noise effects. If the different parameters of 
the algorithm are chosen properly, then the expectation value 
of the energy decreases stably during the updating process. 
Notice that the particular value of the parameters used in the 
procedure also affect the speed and precision of the conver- 
gence. Let us also mention that the VQMC algorithm has a 
more favorable scaling with the MPS parameters than the di- 
rect variational algorithm. 

Once the MPS for the ground state is obtained using any of 
the above two methods, the GE can be computed according 
to its definition. For a system of n spins 1 /2, we maximize 
Eq. d2J with respect to the separable states \<p) of n spins, 
\<P) = ®"=i( cos (£)IO> + sin(£)|l», where |0> and |1) are the 
eigenstates of the S , spin operator with eigenvalue +1/2 re- 
spectively. For the translationally invariant (TI) Hamiltonians 
that we consider in this work, we assume that the maximiza- 
tion of Eq. (O can be done with respect to the set of separable 
states that are also TI, that is, we assume that s £ V; 11911 . 
Let us further justify this assumption: the Hamiltonians that 
we consider in this paper favor a ferromagnetic alignment of 
the spins in the quantum system, that is, spins tend to be par- 
allel to each other. This simplifies the calculation of the geo- 
metric entanglement in that it is possible to optimize over the 
family of TI product states with a period of one site. Notice 
that if the Hamiltonans were such that the spins would tend 
to align antiferromagnetically, then the geometric entangle- 
ment should be obtained by optimizing over a broader family 
of states (most probably over those product states that are TI 
with a period of two sites). As we shall see, this is not the 
case in our work. In fact, for systems of small size, and for the 
Hamiltonians that we will consider here, we have numerically 
checked that the GE obtained by optimizing over the family 
of TI product states with period of one site is equal to that ob- 
tained by optimizing over the family of TI product states with 
period of two sites (see the Appendix), which corroborates the 
validity of our assumption. 

As such, the fidelity A(£) between the separable state \<f>) 
and the MPS approximation to the ground state \tfr g ) takes the 
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form 

A(^) = |<(A g |0>| = |Tr(r [I] r [2] - - ■ r [n] ) |, (8) 

where = £ s .Ag\si) <8> -S(s,) are zero-dimensional transfer 
matrices V/ of some ID classical vertex model, with B(0) = 
cos(£),fi(l) = sin(£), and A g (si) the matrices at site i in the 
MPS representation of the ground state wave function \if/ g ). 
Eq. can be further simplified if the matrices A '''(«,■) are 
site-independent. In this case we have that = T V/ and 
therefore the fidelity can be expressed as A(£) = |Tr(7"")| — 

| 2 fl (/t a )"|, where A a are the eigenvalues of transfer matrix T. 
Finally, the maximization of the fidelity A(f) with respect to 
the angle £ can be achieved by means of standard numerical 
optimization algorithms, and from here the GE per spin £„ 
readily follows from Eqs. (0) and (|4j. Importantly, £„ can be 
understood as twice the free energy per site of the ID classical 
vertex model with partition function A(^ max ), where % ma x is the 
maximizing angle. 

Models. We focus on the quantum Ising spin chain with 
transverse magnetic field, and on the XXZ quantum spin 
chain. The Hamiltonian which accommodates both models 
as two limiting cases takes the form 

H = ~Y i {^-Sfsf I] + ^-Sfsf 1] + jSfSf 1] + /IS™] , 

(9) 

where S„ (en = x, y, z) are the spin operators of the i-th spin 

1 /2, and parameters A, y and A characterize different aspects 
of the model (anisotropies, magnetic field). The Hamiltonian 
from Eq. (O exhibits very rich physics. Here, we consider the 
following two cases: 

(i) For y=l and A = 0, the model corresponds to the 
ID quantum Ising model in a transverse magnetic field, with 
a critical point at A = +1 in the universality class of a free 
Majorana fermionic field theory. Without loss of generality, 
here we consider the case A — 1 . 

(ii) For 7 = 0, A = 0, and A generic, the model is the 
anisotropic XXZ model, which is critical for A € [—1,1] 
and in the universality class of a free compactified bosonic 
field theory. The system exhibits ferromagnetic and anti- 
ferromagnetic orders respectively for A > 1 and A < - 1 l2(ill . 
In this work we focus on the critical regime A e [-1, 1]. 

Simulation results. As a first test, we have compared results 
for the GE obtained using both the direct variational algorithm 
and the VQMC algorithm. In practice, we observe that the di- 
rect variational algorithm seems more stable but breaks trans- 
lational invariance, while the VQMC algorithm imposes trans- 
lational invariance from the beginning but needs more efforts 
to get stabilized. In order to check quantitatively both meth- 
ods, we have computed the GE for the periodic quantum Ising 
chain at the critical point A - I for sizes n from 20 to 28. 
Fixing the precision of the ground state energy in both meth- 
ods to 10" 6 , we see in Table U that the relative error of the 
GE as computed by the two algorithms is always smaller than 

2 x 10~ 5 . Thus, we conclude that both methods are equivalent 



size n 


20 


22 


24 


26 


28 


direct £„ x 10- 


7.7469 


7.3564 


7.0259 


6.7441 


6.4878 


VQMC £„ x 10 z 


7.7464 


7.3563 


7.0263 


6.7436 


6.4884 



TABLE I: The geometric entanglement per spin £„ is calculated by 
means of both the direct optimization and the VQMC algorithms 
for the quantum Ising chain in a transverse magnetic field, at criti- 
cal magnetic field A = 1 and system sizes n from 20 to 28. Fixing 
the precision of the ground state energy at 10 -6 , the results obtained 
from both algorithms for £„ are comparable, with relative differences 
smaller than 2 x 10~ 5 . 
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FIG. 1: (color online) Main: Relation between the geometric entan- 
glement per spin £„ and the chain size n for the ID quantum trans- 
verse Ising model with PBCs. For sizes n from 8 to 120, the data 
are fitted to £„ = a + b/n + f/n 2 . Inset: The relative fitting errors 
| g flt _ £ mpS| /£ mps are always sma n er trian 2.5 x 10~ 3 , where £*' is 
the value extracted from the fit and £j^ ps is the value extracted from 
the MPS algorithm for each n. 

in the calculation of the GE per spin up to a small numerical 
error. 

In Fig. Q]we plot £„ for the quantum Ising model in a trans- 
verse field at the critical point A = 1 as a function of n, for 
sizes from n — 8 up to n — 120. Our data suggest that there is 
indeed a well defined finite-size scaling law behind the behav- 
ior of £„. This finite-size scaling for the GE per spin seems 
to be well fitted by £„ = a + b/n + f/n 2 , where the fitting 
coefficients are a = 0.03096, b = 1.016, and / = -1.713, 
all with relative fitting errors smaller than 2.5 x 10~ 3 . As a 
check, we have also derived the thermodynamic limit of £„ 
from the exact solution of the transverse Ising model in ID 
112 111 , and obtained that at the critical point A = 1 it takes the 
value £oo = 0.03113. This is comparable with the value ob- 
tained from the extrapolation «—> 00 of our fit for £„, which 
is a = 0.03096. Notably, the relative error between the exact 
value £00 and our fitted value a is of the order of 10~ 3 . 

A remark is now in order regarding potential sources of er- 
ror in our fits. First, as n becomes larger, it is clear that the 
fits should be better, therefore large system sizes are required 
in order to get accurate results. Second, the effects of our 
numerical approximation to the ground state based on ma- 
trix product states are also important. Namely, at criticality 
(as is the case) the introduction of a finite bond dimension 
in the matrix product state introduces an artificial correlation 
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FIG. 2: (color online) Main: Relation between the geometric entan- 
glement per spin, £„, and the chain size n for the ID XXZ model 
with PBCs for sizes n from 10 to 94. The plot is for nine different 
values of A > 0. The data are fitted to £„ = a + bjn + f/n 2 . Inset: 
The relative fitting errors |£® - £j) lps |/£jj lps are always smaller than 
lCT 3 , where fijj' is the value extracted from the fit and £j^ ps is the 
value extracted from the MPS algorithm for each n. 



A 


a 


b 


/ 


0.0 


0.15930 


0.9806 


-0.61511 


0.1 


0.14922 


1.0308 


-0.68507 


0.2 


0.13902 


1.0842 


-0.78287 


0.3 


0.12858 


1.1397 


-0.86608 


0.4 


0.11763 


1.2061 


-0.9994 


0.5 


0.10612 


1.2821 


-1.1607 


0.6 


0.09383 


1.3706 


-1.3544 


0.7 


0.08027 


1.4827 


-1.6288 


0.8 


0.06469 


1.6379 


-2.0860 



TABLE II: Fitting coefficients of the geometric entanglement per 
spin, £„, for nine different values of the anisotropy A > in the 
XXZ model. 



length for large system sizes (see the analysis in Ref. [22]). 
This, in turn, moves the ground state of the system slightly 
away from the true critical ground state, which implies that 
our results about critical properties of ID systems can only be 
approximate even if we were able to get the thermodynamic 
limit. The fact that critical systems can only be approximately 
reproduced is a well-known property of all numerical algo- 
rithms based on matrix product states. Our calculations use 
a bond dimension of the matrix product state which is suffi- 
ciently large to assure the correctness of the geometric entan- 
glement. 

We have also computed £„ for the XXZ model with PBCs. 
In the critical regime A e [-1, 1], we choose nine different 
values > of the anisotropy A to investigate the finite size 
behavior of the GE per spin. For sizes n ranging from 10 to 
94 spins, &„ is again well fitted numerically by S„ = a + bjn + 
f/n 2 as seen in Fig. [2] In this case, we obtain different fitting 
coefficients as a function of A, and the relative fitting errors 
are always smaller than 1CF 3 . The specific fitting coefficients 
for each A are given in Table [Hi 

An interesting issue raised by the above results is the pos- 
sible universality of coefficient b in Eq.©. For the quantum 



Ising model, a comparison of our results with those that can 
be extracted from Ref. II 1211 indicates that, indeed, b may actu- 
ally be equal to 1 for this system. For the XXZ model, the fact 
that different values of A correspond to different fitting coef- 
ficients could be an indication that the finite size corrections 
to the GE per site at criticality depend only on the Luttinger 
liquid parameter K of the corresponding effective field the- 
ory at low energies, which in turn labels the corresponding 
universality class. Interestingly, if b were universal then the 
finite-size behavior of the GE would be in sharp contrast to 
that of other quantities such as the ground state energy per 
site, where the universal finite-size scaling correction appears 
at order (9(1 /n 2 ) instead of (9(1 In) II 1 511 - Notice also that other 
entanglement measures such as the von Neumann entropy y|, 
the single copy entanglement and the GE per region for 
a partition into regions of a macroscopic size [5] already dis- 
play universal behaviors. Despite of its inherent interest, the 
findings of this paper are not conclusive about the (potential) 
universality of coefficient b. This problem will be specifically 
addressed in a separate publication 112311 . 

Conclusions and discussion. In this paper we have numeri- 
cally investigated the finite-size GE for ID quantum spin lat- 
tice systems with PBCs in the context of newly-developed TN 
algorithms. In particular, we have shown how to compute the 
GE using these techniques. We have also shown that the lead- 
ing term in the finite-size correction to the GE per spin at crit- 
icality behaves as b/n, as evidenced by our calculations for 
the quantum Ising spin chain in a transverse magnetic field 
and for the XXZ spin chain. The possibility of this leading 
finite-size scaling being universal has also been discussed. 

Although we have restricted ourselves in this work to the 
study of finite size quantum lattice systems in ID, it is pos- 
sible to abandon this restriction in several ways. Specifically, 
one may directly compute the GE per site for finite and infi- 
nite quantum latice systems both in ID and 2D by using TN 
algorithms This will be the subject of future works. 
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APPENDIX 

In this appendix we show a numerical calculation which 
gives evidence that the GE obtained by optimizing over the 
family of TI product states with period of one site is equal 
to that obtained by optimizing over the family of TI product 
states with period of two sites. More specifically, for the XXZ 
model that we consider in this paper with n = 10 sites a simple 
exact diagonalization of the Hamiltonian gives us the results 
from Table (ED. 

Remarkably, in Table ( fTTTb we can see that in both cases (1- 
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A 


GF 1-site TT 


GF 2- site TT 

V 1 1 ' — l31LI^ 1.1. 


0.0 


2.51221 


2.51221 


0.1 


9 45308 


2 45308 


0.2 


2.39471 


2.39509 


0.3 


2 33639 


2 33639 


0.4 


2.28193 


2.28184 


0.5 


2.22632 


2.22632 


0.6 


2.16945 


2.16945 


0.7 


2.11964 


2.12233 


0.8 


2.07288 


2.07288 



TABLE III: GE for the XXZ model with n = 10 sites. The results 
are for 1-site TI product states, and 2-site TI product states. The 
ground state has been obtained by an exact diagonalization of the 
Hamiltonian. 

site TI and 2-site TI) we obtain very similar results for the GE. 
This is an indication that, apparently, the closest product state 
can be chosen to be translationally invariant with period of one 
site, in turn being consistent with our qualitative arguments in 
the main body of the paper. Notice that the above numbers 
also match with the results in Fig. [2] for « = 10 using matrix 
product states. 
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